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Abstract 

A numerical study was performed to investi- 
gate the accuracy and validity of cold-gas simula- 
tion of actual hot scramjet exhaust within a Mach 
10 free stream over a representative si ngle-stage- 
to-orbit airbreathing configuration. In particu- 
lar, exhausts of various noncombusting chemis- 
try models were studied to characterize their 
effects on the vehicle aftbody performance and 
the plume flow field definition. Two approxima- 
tions of the hot scramjet combustion products 
were utilized to determine the requirement for 
expensive, multi-species numerical modeling, 
and to establish a baseline for the validation of 
cold-gas simulation. Cold-gas simulation at 
Mach 10 is shown to be a viable technique using 
an appropriate thermally perfect gas mixture for 
reproducing hot scramjet exhaust effects. 

Introduction 

Prior to the National Aero-Space Plane 
(NASP) program in the United States of America 
very little was known regarding the effects that 
powered scramjet exhaust flow fields would have 
on the aftbody of the vehicle, particularly on the 
aerodynamic and control surfaces. Such exhaust 
flow field effects are commonly referred to as jet- 
effects. Actual air/hydrogen scramjet propulsion 
is not a scalable process, and the si zing of a flow- 
through, fuel-injection, combustion engine sys- 
tem in a typically small hypersonic wind tunnel 
model is extremely difficult. Therefore, other 
approaches, known as exhaust simulation, were 


* Staff Engineer, Senior Member AIAA 
f Aerospace Engineer, Hypersonic Airbreathing Propul- 
sion Branch, Gas Dynamics Division, Senior Member 
AIAA 

This paper is declared a work of the U.S. Government and 
is not subject to copyright protection in the United States 


assessed. Through the course of the NASP pro- 
gram many exhaust simulation techniques were 
studied for obtaining data describing these jet- 
effects using wind tunnel models that produce a 
simulated scramjet exhaust. 

One such approach proposed using a minia- 
ture combustor fed by two separate gas I i nes, one 
for hydrogen and one for oxygen. The resulting 
combustion products approached appropriate 
scramjet exhaust temperatures, but lacked the 
significant amount of nitrogen species found in 
air/hydrogen combustion. Furthermore, this 
technique is very expensive in terms of model 
design and fabrication, as well as in the potential 
expense of modifying a facility for safe operation 
with this type of combustion. 

Alternatively, Oman, et al., 1 developed a 
method to use a non-combusting, thermally per- 
fect gas at 'cold' temperatures tosimulatethe gas 
dynamic effects caused by air/hydrogen combus- 
tion products in a hot exhaust external nozzle 
flow field. Herein 'cold' refers to temperatures 
much less than those occurring in air/hydrogen 
combustion. This simulant (cold) gas concept is 
based upon matching the ratio of specific heats y 
of the exhaust gas, the Mach number, and the 
static pressure ratiothroughout the exhaust flow 
field. Oman, et al., used a mixture of Freon-12 
and argon (Ar), two inert gases, to provide a y 
close to that obtained from actual combustion, 
where y is a function of temperature. Aftbody 
surface pressures obtained using the Freon/Ar 
mixture compared well with combustion data 
obtained on the same model 2 . Because of the 
reduction in facility operations safety concerns 
and model complexity, compared to the air/hydro- 
gen combustion method, this method was 
adopted by NASA Langley for work on studies 
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related to the N ASP 2 3 . (Tetrafluoromethane, CF 4 , 
has since been substituted for the chlorofluoro- 
carbon Freon-12 which has been found to be 
destructive to the ozone layer. CF 4 has proper- 
ties similar to Freon-12, but without the harmful 
chi ori ne component.) 

At the NASA Langley Research Center, the 
Test Technique Demonstrator (TTD) model was 
conceived to investigate a wide variety of wind 
tunnel testing techniques associated with air- 
breathing, single-stage-to-orbit (SSTO) vehicles. 
As a complementary effort, computational fluid 
dynamics (CFD) solutions of this configuration 
were performed for CFD code calibration/valida- 
tion and to better understand the flow physics 
associated with this class of vehicle. A large part 
of this experimental and computational research 
was focussed on understanding the effects of 
scramjet exhaust aftbody flow fields. Computa- 
tionally, these included powered effects due to 
angle of attack 4 , and effects of inlet representa- 
tion on aftbody flows 5 . I n the current study, dif- 
ferences in aftbody performance and plume flow 
field definition resulting from various degrees of 
complexity of exhaust gas modeling were exam- 
ined computationally through the use of four 
exhaust chemistry models. Previous studies at 
Mach 6 6,7 concluded that cold-gas simulation of 
scramjet exhaust flows is a viable approximation 
to actual hot scramjet exhaust at that Mach num- 
ber. The cold, reduced-y exhaust of a CF 4 -Ar mix- 
ture yields slightly better simulation of the 
pi ume size and shape, and the aftbody forces and 
moments, due to scramjet exhaust than does a 
cold air exhaust at Mach 6. 

Flowever, the preceding studies did not 
address the question of the Mach number range 
for which the cold-gas simulation technique is 
valid. I n particular, does the technique retain its 
accuracy relative to hot exhaust flows in the free 
stream Mach number range (approximately 6-15) 
in which scramjets are primarily intended to 
operate? The current paper describes a numeri- 
cal study of various cold and hot exhaust simula- 
tions at Mach 10. The objectives of this study 
were to determine the accuracies of aftbody/ 
pi u me effects due to col d-gas si mu I at i ons rel at i ve 
to hot (noncombusting) scramjet exhaust, includ- 
ing two approximations to the hot scramjet 
exhaust. 


Modeling Considerations 

Model Geometry 

The Test Technique Demonstrator (TTD) is a 
generic N ASP-1 ike configuration designed prima- 
rily for investigating wind tunnel test techniques 
for hypersonic air-breathing vehicles. Two small 
scale models (approximately 1.5 and 3 feet long) 
have been fabricated for scramjet exhaust simu- 
lation powered testing in NASA Langley’s high- 
speed wind tunnels 8,9 . Reference 8 describes 
methodology for using metric model parts to 
obtain force and moment data under powered 
conditions, and reference 9 describes the acquisi- 
tion of powered effects data using surface and 
flow field pressure measurements. The concep- 
tual airbreathing flight vehicle configuration is 
shown in F igure 1. H owever, as noted previously, 
true powered airbreathing, fuel-injection, com- 
bustion testing is impractical in current small- 
scale high-speed wind tunnels due to both tunnel 
and model size constraints. I nstead, the forebody 
inlet is closed with a fairing and exhaust simula- 
tion is performed. The fairing eliminates the 
problem of removing the air captured by the 
inlet, thereby allowing a simulant gas to be 
routed on board and expanded out the nozzle. 
Computational studies have shown that a prop- 
erly designed fairing has minimal effect on the 
aftbody/exhaust flow field 5,10 . A lower aft view of 
the fai red-over configuration modeled in the cur- 
rent study isshown in Figure2, all owing the two- 
module internal nozzle to be seen. Figure 3 
shows how a simulant exhaust gas is routed on 
board the actual wind tunnel model through the 
model support strut into a plenum chamber and 
exhausted out an internal nozzle designed to 
yield specified flow conditions at the exit, as 
described below. The support strut and external 
gas supply system werenot modeled inthis work. 
The internal nozzle is contained within the 
engine module as part of the forebody, and the 



Figure 1. Test Technique Demonstrator 
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Figure 2. Full-span shaded numerical model of theTTD with wings and inlet fairing 
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Figure 3. Schematic of small-scale model exhaust simulation flowpath 


nozzle exit plane (engine cowl trailing edge) 
defines the downstream end of the forebody. The 
nozzle exit plane is designated station 4, consis- 
tent with previous report terminology. The aft- 
body is the remainder of the vehicle, i.e., 
everything aft of the cowl trailing edge (nozzle 
exit plane) and includes the wings, mounted at a 
nominal incidence of -1.5° (leading edge down) for 
this study. 

The simulated scramjet exhaust in this 
numerical study was computed beginning with a 
uniform Mach number of 1.03 just downstream of 
the nozzle throat and expanded to the cowl trail- 
ing edge (see Figure 3). Due to the variety of 
exhaust gas chemistry models employed in this 
study, each having different flow expansion char- 
acteristics, four different internal nozzles were 
requi red to provide average cowl exit-plane M ach 
numbers of approxi mately 3.4 for the appropri ate 


exhaust gas. The nozzles for hot exhaust gases 
were designed to also yield average exit-plane 
static temperatures of approximately 3410°R 
(1895 K), typical of actual scramjet exit tempera- 
tures. The exit plane geometry was held fixed 
and the nozzle length was allowed to vary only 
0.2 inches. The nozzle walls had sharp trailing 
edges to avoid recirculation regions aft of the 
cowl, thus allowing supersonic space-marching 
CFD solutions. An internal sharp-edged splitter 
plate at the vehicle symmetry plane divided the 
nozzle into two modules. Figure 4 shows a three- 
view schematic of a typical internal nozzle mod- 
ule. 

Flow Conditions 

The purpose of this study was to analyze the 
validity of cold-gas powered simulation at a typi- 
cal flight Mach number in which a scramjet 
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Figure 4. Internal nozzle schematic 


engine would be operating. A baseline nozzle 
pressure ratio (NPR = Pt/pJ of 3000 was 
selected as representative of Mach 10 flight, 
where the 't,j ' refers to total (plenum) conditions 
of the exhaust jet. The free stream conditions 
were selected to match nominal conditions in the 
NASA LaRC 31-Inch Mach 10 Wind Tunnel: 
Re^ = 2xl0 6 /ft, = 10.14, p t oo = 1450 psi, and 

T t|CO = 1836°R. I n all cases, the external gas was 
thermally perfect air, modeled as a single species. 
A single forebody CFD solution based upon these 
free stream conditions was computed and com- 
bined with various internal nozzle solutions to 
initializethe aftbody flow field for each solution. 

Twelve internal nozzle/exhaust combinations 
were computed using four noncombusting, ther- 
mally perfect gas mixtures at both cold and hot 
(combustion) temperatures (see Table 1). Four 
cases used air as the exhaust gas, three at cold 
conditions for a range of NPR values and one hot 
case at the baseline NPR. Three cases used a 70/ 
30% mole fraction mixture of CF 4 and Ar at cold 
temperatures, the current standard simulant gas 
in use at NASA LaRC. These exhaust calcula- 
tions were compared with two approximations to 
the actual composition of Mach 10 scramjet com- 
bustion products at hot temperatures. A six-spe- 
cies approximation, consisting of about 30% H 2 0 
(steam), 63% N 2 , 1% 0 2 , 1% NO, 1% Ar, and 4% 
H 2 , by molefraction, accounts for all but residual 
species from the actual combustion process. A 
two-species model of approximately 32% H 2 0 
and 68% N 2 by mole fraction preserves the rela- 
tive concentrations of the two dominant combus- 


Table 1: CFD Case Summary 


Exhaust gas 

NPR 

T t :j (°R) 

Air 

2050 

790 

Air 

2700 

790 

Air 

3000 

790 

Air 

3000 

9360 
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3000 

790 
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3300 

790 

L_ 
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3600 

790 

h 2 o-n 2 

2700 

8460 

h 2 o-n 2 

3000 

8460 

h 2 o-n 2 

3300 

8460 

H 2 0-N 2 -0 2 -H 2 -N0-Ar a 

3000 

8460 

H 2 0-N 2 -0 2 -H 2 -N0-Ar a 

3300 

8460 


a. same nozzle as H 2 0-N 2 . 


tion product species at a lower computational 
cost. Five calculations were performed with 
these approximations to determine NPR effects 
for the hot exhaust and to determi ne the requi re- 
mentsfor precisely modeling the mixture of com- 
bustion products. The same nozzle geometry and 
jet total temperature was used for all of the hot 
combustion products exhaust calculations. 

While the NPR is an appropriate figure of 
merit for scramjet engine designers, the impor- 
tant simulation parameters 1 are the local condi- 
tions at the nozzle exit. For more appropriate 
comparisons of the resulting data, a static nozzle 
pressure ratio (SNPR 4 = p^pj was determined 
from the computed flow fields, where p 4 was eval- 
uated just upstream of the cowl exit plane on the 
body-side surface of the internal nozzle (the loca- 
tion of pressure taps in the actual wind tunnel 
models). 

All hot calculations were made under the 
assumption of thermally perfect, frozen flow, i.e., 
no combustion or dissociation. This assumption 
was reasonable downstream of the nozzle/cowl 
exit plane wherescramjet combustion is typical ly 
considered to be complete and the static temper- 
atures were below 3500°R, even though it was 
obviously not valid at the extremely high internal 
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nozzle temperatures. Since the external flow 
fields were the emphasis of this study, the (ficti- 
tious) internal flows were computed only to pro- 
vide reasonable initial exhaust flow fields to the 
aftbody. The six and two-species combustion 
products compositions were assumed to be those 
at the exit plane, as provided by the frozen inter- 
nal flows. 

Each of the internal nozzle solutions was 
merged with the single forebody solution via 
bilinear interpolation at the cowl exit plane to 
initialize each aftbody flow field solution. The 
cases comprisi ng the study (Table 1) were charac- 
terized by the exhaust gas mixture, jet condi- 
tions, and a corresponding internal nozzle 
geometry. Figure 5 shows the relationships 
between the NPR and SNPR 4 for each case. 
(Note that the NPR values for the hot cases are 
valid only for computational purposes under the 
assumption of frozen flow. Experimental NPR 
values would be different due to the chemical 
reactions and dissociation occurring within the 
internal nozzle.) The two hot combustion prod- 
ucts exhausts agree to within 0.6%. At the base- 
line NPR of 3000, the cold air exhaust is 43% 
higher in SNPR 4 compared to the most complete 
(six-species) combustion products model, while 
the cold CF 4 -Ar exhaust is 4.7% lower. The 
slopes of the SN PR 4 versus NPR for the two hot 
combustion products exhausts agree within 
about 1.5%. The cold air exhaust slope is 41% 
high, while the CF 4 -Ar simulant slope is 8% low. 


This wide range of data slopes complicates the 
task of choosing simulant NPR values to achieve 
appropriate station 4 conditions, especially for 
cold air. Using hot air does reducethe magnitude 
of the difference in SNPR 4 considerably, but this 
isnomoreviableasatesttechniquethan is using 
hot combustion products. The hot air exhaust 
data is presented only to illustrate temperature 
effects while holding the gas fixed. 

CFD Code Description 

The General Aerodynamic Simulation Pro- 
gram (GASP), version 2 . 0 U , was the CFD code 
applied in the current study. GASP solves the 
integral form of the governing equations, includ- 
ing the full Reynolds-averaged Navier-Stokes 
equations, theThin-Layer Navier-Stokes (TLNS) 
equations, the Parabolized Navier-Stokes (PNS) 
equations, and the Euler equations. A general- 
ized chemistry model and both equilibrium and 
non-equilibrium thermodynamics models are 
included. The discretized equations may be 
solved by space marching or global iteration, and 
both explicit multi-stage Runge-Kutta and 
implicit time integration schemes are included. 
Characteristic-based spatial discretizations 
employ a variety of flux and J acobian evaluation 
schemes. 

The code allows highly flexible coupling of 
multiple grid zones, each of which may be solved 
independently or in a coupled manner, with spe- 
cies differences permitted between zones, as well 



Figure 5. Variation of static nozzle pressure ratio with total nozzle pressure ratio 
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as choices of numerical algorithms. All thermo- 
dynamic and transport property data for individ- 
ual species are permanent parts of the code, 
selected at run time through input options. The 
solution algorithms, flux models, limiters, vis- 
cous models, boundary conditions, and initial 
flow conditions are all set for each zone within 
the input files. Grid generation for GASP is sim- 
plified due tothe zonal algorithm, but grid points 
of adjacent zones must coincide along grid inter- 
face boundaries for global, or single-run, compu- 
tations. 

Before any CFD code can be used as an analy- 
sis tool, code calibration is essential for the types 
of flows of i nterest to provide a level of confidence 
in the ability of the code to accurately predict the 
fluid dynamics of the flow fields. Previous studies 
using GASP have shown that it has the ability to 
accurately predict complex, three-dimensional 
hypersonic flows including a representative 
NASP forebody 12 , the effects of plume-induced 
separation 13 , and three-dimensional powered 
effects on a model aftbody 5 . 

Numerical Issues 

All CFD calculations in this study solved the 
PNS equations using third-order upwind-biased 
inviscid split fluxes. Laminar flow was assumed 
based upon the small model scale and the rela- 
tively low free stream Reynolds number. The 
thermodynamics model made use of the NASA 
Lewis curve-fit 14,15 forms for thermally perfect 
gases, with local extensions of the valid ranges 
derived from J ANAF 16 and National Bureau of 
Standards 17 data. Such extensions were neces- 
sary to model the I ow free stream static tempera- 
tures representative of the Langley 31-Inch 
Mach lOTunnel. The transport properties were 
based upon Sutherland's model, using Wilke's 
rule for mixtures. 

Wall temperatures on all model surfaces were 
held constant, as if the model was being actively 
cooled. For the majority of cases thewall temper- 
ature^ was chosen as 90°F, but the GASP code 
exhibited numerical instability for the two cases 
employing the hot six species exhaust gas model. 
For these cases the wall temperature was 
increased to 260°F, the minimum to achieve 
numerical stability. The effect of this change in 
T w remai ns to be determi ned. 

Separate structured grids were generated for 
theTTD forebody, each internal nozzle, and the 


aftbody. Only half of the full configuration was 
modeled based upon symmetry about a vertical 
center-pl ane. The forebody grid di mensi ons were 
45 cells streamwise, 64 cells normal tothe body, 
and 64 cells wrapping around the body from 
upper symmetry plane to lower symmetry plane. 
The internal nozzle grids all consisted of 44 cells 
streamwise, 80 cells vertically, and 68 cells span- 
wise from the center splitter pi ate to the si dewall. 
The aftbody grid dimensions were 70 cells 
streamwise, 64 cells normal tothe body, and 128 
cells wrapping from upper to lower symmetry 
plane. Thus, the grid sizes for the half-span 
numerical models were forebody: 184,320 cells, 
nozzles: 239,360 cells, and aftbody: 573,440 cells, 
for a sum total of 997,120 computational cells. 
The forebody calculation took about 775 cpu sec- 
onds on a CRAY C-90. Typical internal nozzle 
solutions required 300 C-90 cpu seconds for cold 
air and 2600 seconds for a hot six-species 
exhaust. Typical aftbody calculations required 
1100 C-90 cpu seconds for cold air exhaust, 3700 
seconds for hot two-species exhaust, and 21,550 
seconds for hot six-species exhaust. 

The grids were clustered near the body for 
accurate resolution of boundary layers. The 
inner law variable, y + = Y [ /p 1 , provides 

a measure of accuracy of the viscous sol utions by 
combining density p, velocity u, and viscosity p 
with the amount of point spacing normal to the 
surface AY, all computed at the first computa- 
tional cell center off the body (denoted by the sub- 
script 1). Typically, y + values of the order of lwill 
ensure reasonable grid- point resolution of turbu- 
lent shear layers and, thus, excellent resolution 
of laminar boundary layers. I n the forebody cal- 
culation the cell y + values averaged 0.76, with a 
local maximum of 4.3, occurring along the out- 
board edge of the inlet fairing. On most of the 
fairing surface they + values averaged about 2. I n 
the aftbody calculations the average y + values 
were between 2 and 2.5, with local peaks of 6.7 at 
the lower wing/body juncture and a spike of 24 
along the wing leading edge (near the stagnation 
line). 

External body forces and moments were cal- 
culated by integration of computed surface pres- 
sure and skin friction data. Bilinear 
distributions of these quantities were assumed 
over each quadrilateral surface grid cell. Normal 
force components were computed by integration 
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of the pressure over projected areas normal to 
each of the three cartesian directions. Tangential 
force components were computed from surface 
skin friction coefficients using velocity compo- 
nents at the first cell center off the body surface 
to resolve directions. The normal and tangential 
force components were summed over all surface 
grid cells to obtain total force components. Cell 
pitching moments were computed assuming that 
the cell forces acted at the cell centroid. Summa- 
tion of the cell pitching moment values yielded 
total body pitching moments. 

Results 

Data comparisons presented below consist of 
two basic types: external forces and moments on 
the aftbody, and contours of plume parameters, 
in order to illustrate both surface and flow field 
effects. 

Forces and Moments 

Figure 6 shows longitudinal aftbody force and 
moment components normalized by the magni- 
tudes of the forces and moments of the hot scram- 
jet combustion products simulation, using the 
(more complete) six-species combustion products 
mixture at the baseline NPR=3000, plotted ver- 
sus SNPR 4 . Each of the simulant exhaust 
approximations described in Table 1 are shown, 
including the single data point for hot air 
exhaust. I n general, for a given exhaust mixture 
at a given total temperature the variations with 
respect to SNPR 4 are nearly linear. Note the 
I argeSNPR 4 value (50.8) for the cold air exhaust 
at the baseline NPR=3000 compared to the other 
gases (see Figure 5). The NPR=2050 cold air case 
was run to allow comparisons at the SNPR 4 of 
the baseline case (SNPR 4 =35.5) without extrapo- 
lation. 

The aftbody lift shows the least sensitivity to 
exhaust mixture or temperature. The slopes of 
the curves all agree to within 2%. The two-spe- 
cies approximation of combustion products 
agrees well with the six-species model. The cold 
air results show excellent agreement at corre- 
sponding values of SNPR 4 , whilethecold CF 4 -Ar 
results are 4.8% I ow. The pressure forces account 
for greater than 99% of the total aftbody lift 
regardless of exhaust gas; i.e., the contribution of 
skin friction tothe lift is negligible. 

Aftbody thrust is noticeably more sensitive to 
gas composition and temperature. Similar to lift, 
the slopes of thrust with respect to SNPR 4 agree 
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Figure 6. TTD aftbody force and moment 
comparisons 

to within 2% except for cold air, which is approx- 
imately 10% greater than that of the six-species 
combustion model. The cold air thrust data over 
predicts the six-species combustion products 
thrust by about 21% at the baseline SNPR 4 . The 
cold CF 4 -Ar thrust data are about 6% higher, and 
the two-species combustion model 6% lower, than 
the six-species combustion model. Theskin fric- 
tion contribution (actually, drag) amounts to a 
significant percentage of total aftbody thrust. 
Specifically, the viscous drag is about 27% of the 
magnitude of the pressure contribution for the 
hot six-species exhaust, and about 20% for the 
cold CF 4 -Ar exhaust. This decrease in viscous 
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drag, relative to the pressure thrust, for the cold 
exhausts contri butes to the correspond! ng higher 
aftbody thrusts. 

The trends in pitching moment are combina- 
tions of the trends in lift and thrust since the 
moment calculations are functions of both force 
components. The combustion models agree in 
magnitude within about 1.5%, and the cold CF 4 - 
Ar and air results are 8% and 6% high, respec- 
tively. The slopes of pitching moment versus 
SNPR 4 all agree within about 4%. Similar to lift, 
the pressure contri bution is about 99% of the aft- 
body pitching moment. 

The single hot-air data point on each plot of 
Figure 6 illustrates the temperature effects on a 
specific gas. The combustion products mixtures 
and the cold CF 4 -Ar mixture have ratios of spe- 
cific heats in the range of 1.25 to 1.28 at the noz- 
zle exit. At nozzle exit temperatures of 3400°R 
the ratio of specific heats for air is reduced from 
1.4 at 'cold' temperatures to approximately 1.3. 
Flowever, within the internal nozzle the hot-air 
temperatures are even higher, dissociation would 
be expected, and the thermally perfect assump- 
tion for air made in this ideal computational 
study is no longer valid. Thus, heating air 
exhaust to high temperatures to reduce the ratio 


of specific heats is not practical, even though the 
ideal thermally perfect data shown herein would 
indicatethat results from such a technique would 
be better than using cold air. 

Plume Flow Field 

Figure 7 illustrates the resultant exhaust 
plume shapes and extents for two of the solutions 
at similar values of SNPR 4 . J et exhaust mass 
fraction contours are plotted on a plane normal to 
the free stream at the body trailing edge. The 
right half of the plot isthe hot six-species exhaust 
model of the scramjet combustion products. The 
left side shows similar contours for a cold-gas 
simulation using the CF 4 -Ar exhaust mixture. 
The cases are characterized by the exhaust gas 
mixture and by the value of T 4 , an average tem- 
perature across the nozzle exit plane (station 4). 

I n general, the plumes aresimilar in both size 
and shape. The cold CF 4 -Ar plume is slightly 
larger both vertically and laterally, but at a 
slightly higher value of SNPR 4 (about 3.3%). 
Thus, some of the plume size difference is proba- 
bly due to the increase in SNPR 4 . Figure 6 shows 
a corresponding 5.7% increase in thrust for the 
cold CF 4 -Ar exhaust. Flowever, the same figure 
shows a 4.3% decrease in lift for the cold CF 4 -Ar 



Figure 7. Plume exhaust mass fractions at the TTD body trailing edge 
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exhaust, indicating that pressure gradients 
within the plume are affected by the exhaust 
mixture. Figure 8 supports these trends and 
illustrates this change in pressure gradients by 
comparing pressure ratio contours (p/pj on the 
lower surface of the aftbody for the cold CF 4 -Ar 
and the hot six-species exhausts. Note along the 
centerline the reversal in the mismatch of con- 
tours 2 and 10. Figure 8 also shows that, for 
either exhaust gas, the entire lower aftbody sur- 
face has a positive net force (i.e., no overexpan- 
sion). A similar analysis of the skin friction (not 
presented) shows that viscous drag of the cold 
exhausts is about half that of the hot exhaust. 

The simulant gas concept 1 is ideally based 
upon matching Mach number, static pressure 
ratio, and y throughout the exhaust flow field. 
These three parameters were closely achieved 
within this study, except for the y of the air 
exhaust. All of the internal nozzles were 
designed for an exit plane Mach number of 
approximately 3.4, and SNPR 4 values of approx- 
imately 35 were achieved for all exhaust gases. 
Figure 9 shows the variation of y along an arbi- 
trary streamwise computational grid line within 


the plume for five calculations plotted versus a 
normalized aftbody length (cowl trailing edge to 
body trailing edge). As previously noted, they of 
cold air is a poor match for that of the combustion 
products. The hot air y is somewhat better, but 
still high. The cold CF 4 -Ar agrees much better at 
the nozzle exit, but the gradients of y versus 
length are too large just downstream of the cowl 
trailingedge. This variation issimilartothesur- 
face pressure gradient variation seen in Figure8. 

Thus, the high pressure exhaust flow at the 
nozzle exit expands more rapidly in the CF 4 -Ar 
plume than in the combustion products plume 
due to the thermodynamic properties of the gas. 
The overall plume size and shape is affected only 
slightly, as is the thrust which is dominated by 
the pressures near the nozzle. FI owever, the aft- 
body lift force is distributed more evenly over the 
entire aftbody and is affected more by the 
streamwise pressure gradients. 

Conclusions 

A computational study has been performed to 
investigate the accuracy and validity of cold-gas 
simulation of actual hot scramjet exhausts 


CF 4 -Ar Exhaust 
T 4 = 195°R 
SNPR, = 40.4 


6-species Exhaust 
T 4 = 3360°R 
SNPR, = 39.1 



Figure 8. TTD aftbody lower surface pressure ratio contours, p/p^ 


9 

American Institute of Aeronautics and Astronautics 



AIAA 95-6068 



Normalized Aftbody Length (cowl trailing edge = 0.0) 

Figure 9. Comparison of exhaust gas ratios of specific heats within the exhaust plume 


within a Mach 10 free stream over a representa- 
tive single-stage-to-orbit airbreathing configura- 
tion. In particular, exhausts of various 
noncombusti ng chemistry models were studied to 
characterize their effects on the vehicle aftbody 
performance and plume flow field definition. Two 
approximations of the scramjet combustion prod- 
ucts were utilized to determine the requirement 
for expensive, multi-species numerical modeling, 
and to establish a baseline for the validation of 
cold-gas simulation. 

Based solely upon economy, cold (ambient) air 
is the optimum simulant exhaust gas. It is 
readily avail able at high pressure and poses little 
or no environmental or biological hazards. How- 
ever, from a technical viewpoint, cold air was 
shown to provide poor simulation of the aftbody 
thrust. Hot air appears to show improved thrust 
results relative to the combustion products, but 
the exhaust temperatures required are higher 
than those at which air begins to dissociate. 
Wind tunnel testing at such temperatures would 
be difficult and the air exhaust would not bether- 
mally perfect. 


A cold CF 4 -Ar mixture provides a much better 
representation of the hot combustion products 
exhaust. The ratio of specific heats at the nozzle 
exit more closely matches that of hot scramjet 
exhaust than does they of cold or hot air. The lift, 
thrust, and pitching moment components are all 
simulated within approximately 5-6% accuracy 
by the cold CF 4 -Ar exhaust whiletheair exhaust 
accuracy is inconsistent: good for lift, poor for 
thrust, and moderate for pitching moment. 

The plume definition of the cold CF 4 -Ar 
exhaust agrees well with the combustion prod- 
ucts exhaust. Differences of pressure gradients 
within the plume are shown to be functions of 
thermodynamic properties of the two mixtures. 
H owever, these gradient differences are of much 
less importance in scramjet exhaust simulation 
than the basic premises of matching static condi- 
tions at the nozzle exit, particularly y. It is possi- 
blethat a better simulant exhaust gas than a 70/ 
30% mole fraction mixture of CF 4 and argon can 
be identified to improve the simulation of actual 
scramjet combustion products with regard to y. 
These results show that cold-gas simulation at 
Mach 10 is a viable technique for reproducing 
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hot scramjet exhaust effects when an appropriate 
thermally perfect gas mixture is used. 

Finally, the use of two different hot approxi- 
mations of the scramjet combustion products in 
the numerical solutions showed that the less 
abundant components (4% and less, each, by 
mole fraction) do affect the exhaust characteris- 
tics slightly. However, only the thrust shows a 
discrepancy of more than 2%. Thus, thetwo-spe- 
cies approximation is acceptable if the computa- 
tional cost is significant, but for accurate thrust 
calculations the six-species model should be 
employed. 

References 

1. Oman, R. A.; Foreman, K. M .; Leng, J .; and 
Hopkins, H. B.: Simulation of Hypersonic 
Scramjet Exhaust. NASA CR-2494, March 
1975. 

2. Hopkins, H. B.; Konopka, W.; and Leng, J .: 
Validation of Scramjet Exhaust Simulation 
Technique. NASA CR-2688, J une 1976. 

3. Cubbage, J ames M .; and Monta, William J .: 
Parametric Experimental I nvestigation of a 
Scramjet Nozzle at Mach 6 With Freon and 
Argon or Air Used for Exhaust Simulation. 
NASA TP 3048, February 1991. 

4. Tatum, Kenneth E.; and Huebner, Lawrence 

D. : A Preliminary CFD Study of Angle-of- 
Attack Effects on a Powered Hypersonic 
Model. NASP CR 1154, April 1993. 

5. Huebner, Lawrence D.; and Tatum, Kenneth 

E. : CFD Code Calibration and Inlet-Fairing 
Effects on a 3D Hypersonic Powered-Simula- 
tion Model. AIAA-93-3041, AIAA 24th Fluid 
Dynamics Conference, J uly6-9, 1993. 

6. Tatum, Kenneth E.; Huebner, Lawrence D.; 
and Witte, David W.: Computational Study on 
the Use of Simulant Gas in Powered, Hyper- 
sonic Air-Breathing M odds. NASP TP 1007, 
October 1992. 

7. Tatum, Kenneth E.; and Huebner, Lawrence 
D.: A 3D CFD Study on Simulated Exhaust 
Plumes. Paper No. 38, 1993 National Aero- 
Space Plane Technology Review, Monterey, 
CA, April 13-16, 1993. 

8. Witte, David W.; H uebner, Lawrence D.; and 
Haynes, Davy A.: Test Technique for Using 
MetricModd PartsTo Obtain Powered Effects 


on Air-Breathing Configurations in Hyper- 
sonic Facilities. NASP TP 1008, October 
1993. 

9. Everhart, J oel L.; Ashby, George C., J r.; and 
Monta, WilliamJ .: Test Description and Pre- 
liminary Pitot-Pressure Surveys for Langley 
Test Technique Demonstrator at Mach 6. 
Al AA-92-3940, AIAA 17th Aerospace Ground 
Testing Conference, J uly 6-8, 1992. 

10. Huebner, Lawrence D.; and Tatum, Kenneth 
E.: Computational Effects of Inlet Represen- 
tation on Powered Hypersonic, Airbreathing 
Models. AIAA J ournal of Aircraft, Vol. 30, No. 
5, September-October 1993, pp. 571-577. 

11. McGrory, William F.; Huebner, Lawrence D.; 
Slack, David C.; and Walters, Robert W.: 
Development and Application of GASP 2.0. 
Al AA-92-5067, December 1992. 

12. Huebner, Lawrence D.; and Haynes, Davy 
A.: Forebody Redesign and Flow Character- 
ization of Test Technique Demonstrator at 
Mach 6. NASP TP 1013, April 1994. 

13. Huebner, Lawrence D.; and Tatum, Kenneth 
E.: Computational and Experimental Aftbody 
Flow Fields for Hypersonic, Airbreathing 
Configurations with Scramjet Exhaust Flow 
Simulation. Al AA-91-1709, AIAA 22nd Fluid 
Dynamics, Plasma Dynamics & Lasers Con- 
ference,) une 24-26, 1991. 

14. McBride, Bonnie J .; Heimel, Sheldon; 
Ehlers, J anet D.; and Gordon, Sanford: Ther- 
modynamic Properties to 6000°K for 210 Sub- 
stances Involving the First 18 Elements. 
NASA SP-3001, 1963. 

15. Rate Constant Committee, NASP High- 
Speed Propulsion Technology Team: Hyper- 
sonic Combustion Kinetics. NASP TM-1107, 
NASP J PO, Wright-Patterson AFB, 1990. 

16. J ANAF Thermochemical Tables, Second ed. 
NSRDS-NBS 37, U. S. Dep. Commerce, J une 
1971. 

17. Hilsenrath, J oseph; Beckett, Charles W.; 
Benedict, William S.; Fano, Lilia; Hoge, 
Harold J .; Masi, Joseph F.; Nuttall, Ralph L.; 
Touloukian, Yeram S.; and Woolley, Harold W.: 
Tables of Thermal Properties of Gases. N BS 
Circ. 564, U. S. Dep. Commerce, Nov. 1, 1955. 


11 

American Institute of Aeronautics and Astronautics 



